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■ ABSTRACT 

i ■ 

Q . Context. The majority of galactic chemical evolution models assumes the instantaneous mixing approximation (IMA). This assumption is 
- probably not realistic as indicated by the existence of chemical inhomogeneities, although current chemical evolution models of the Milky Way 
can reproduce the majority of the observational constraints under the IMA. 

Aims. The aim of this paper is to test whether relaxing this approximation in a detailed chemical evolution model can improve or worsen the 
agreement with observations. To do that, we investigated two possible causes for relaxing of the instantaneous mixing: i) the "galactic fountain 
I , time delay effect" and ii) the "metal cooling time delay effect". 

■ Methods. We considered various galactic fountain time delays for the chemical enrichment from massive stars. We then tested a time delay in 
C~~) " the enrichment from stars of all masses due to gas cooling in the range 0.5-1 Gyr. 

f^) , Results. We found that the effect of galactic fountains is negligible if an average time delay of 0.1 Gyr, as suggested in a previous paper, is 

■ assumed. Longer time delays produce differences in the results but they are not realistic. We also found that the O abundance gradient in the 
disk is not affected by galactic fountains. The metal cooling time delays produce strong effects on the evolution of the chemical abundances 

\Q . only if we adopt stellar yields depending on metallicity. If instead, the yields computed for to the solar chemical composition are adopted, 
f^**) ' negligible effects are produced, as in the case of the galactic fountain delay. 

0^ , Conclusions. The relaxation of the IMA by means of the galactic fountain model, where the delay is considered only for massive stars and 

■ only in the disk, does not affect the chemical evolution results. The combination of metal dependent yields and time delay in the chemical 
^ enrichment from all stars starting from the halo phase, instead, produces results at variance with observations. 
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1. Introduction 

The first attempts to follow the chemical evolution of galax- 
ies (Schmidt 1963, Searle & Sargent 1963, Tinsley 1974) were 
based on the Simple Model assumptions: 1) the system is one- 
zone and closed, namely there are no inflows or outflows, 2) the 
initial gas is primordial (no metals), 3) the initial mass function 
(IMF) is constant in time, 4) the gas is well mixed at any time 
(instantaneous mixing approximation or IMA), 5) stars more 
massive than 1 M die instantaneously, while stars smaller than 
1 M live forever (instantaneous recycling approximation or 
IRA) . While the IRA has been relaxed in most of the recent 
chemical evolution models, there have been only a few attempts 
to relax the IMA mainly because of difficulties in estimating 
the dispersal and mixing time of the chemical elements (Roy 
& Kunth 1995) and because models which retain IMA pro- 
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vide nevertheless an excellent fit to the data in the Milky Way 
(e.g. Francois et al. 2004, hereafter F04). Although there is no 
doubt that IMA is not realistic (Schmidt 1963; Tinsley 1975), 
the validity of this assumption depends on the time-scale of 
the mixing processes. Malinie et al.(1993) claimed that due to 
chemical inhomogeneities in the disk, mixing and star forma- 
tion might be delayed by 10 8 ~ 9 yr. 

One of the very few studies in which IMA has been relaxed 
is the one of Thomas et al. (1998). They proposed a model in 
which relaxing IMA was included by means of the splitting of 
the gas components into two different phases: the inactive and 
active gas phases. The first one consists in a hot and not ho- 
mogeneously distributed component; stars cannot form out of 
this phase. The active phase is instead assumed to be cool and 
well mixed. They discussed various time-scales for the mixing 
of the stellar ejecta, and they found that a delay of the order of 
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10 years leads to a better fit of the observational data in the 
[Mg/Fe]-[Fe/H] diagram, than in the IMA hypothesis. 

In this paper we test the relaxation of IMA in the detailed 
chemical evolution model of F04, considering two physically 
motivated time delays: i) the "galactic fountain time delay" ef- 
fect and ii) the delay produced by the cooling of the metals 
("metal cooling delay model"). 

Type II SNe usually occur in OB associations containing 
several dozen massive stars. Sequential supernova explosions 
create supershells which can break out of a stratified medium, 
producing bipolar outflows. The gas of the supershells can frag- 
ment into clouds which eventually fall toward the disk produc- 
ing the so-called galactic fountains. Spitoni et al. (2008) (here- 
after Paper I) followed the clouds ejected by SN explosions in 
the Galaxy both in the purely ballistic model (the cloud is sub- 
ject only to the gravitational potential of the Galaxy) and in 
the hybrid one where viscous interactions between the extra- 
planar halo gas and the cloud are taken into account. They 
found that the range of the cloud orbits is quite small. The 
clouds are generally directed outwards but the average landing 
coordinates differ from the throwing coordinates by at most ~ 1 
kpc. Only for a throwing coordinate of 12 kpc and an OB asso- 
ciation made of 500 SNe the ballistic model predicts a landing 
coordinate ~ 2 kpc larger than the throwing one. This result 
is consistent with the works of Bregman (1980), Fraternali & 
Binney (2008) and Melioli et al. (2008a). Melioli et al. (2008b) 
presented results for a multiple fountain model produced by 
randomly clustered explosions of SNe originating in stellar as- 
sociations spread over a disk area of 8 kpc 2 . As in the case of 
a single fountain, the spreading of the SN ejecta back into the 
disk is not very large. Most of the gas lifted up by fountains 
falls back within a distance AR=±0.5 kpc. Therefore it is un- 
likely that galactic fountains can affect abundance gradients on 
large scales. 

However, as shown in Paper I, the fountains take a finite 
and non-negligible time to orbit around and fall back into the 
Galaxy. This implies a delay in the mixing of metals in the 
interstellar medium (ISM), which conflicts with the IMA as- 
sumption. Referring to the results of Paper I, in this paper we 
test the effect on the chemical evolution of different values of 
the delay originated by sequential SN explosions at different 
galactocentric distances. 

In this work we also investigate another possible cause of 
delayed mixing of the ISM such as the delay due to the fact that 
stars form only in cold gas and that hot gas takes a finite time 
to cool. However, there is still some debate in the literature on 
the cooling time-scale. Recchi et al. (2001) studied the effect of 
a single, instantaneous star-burst on the dynamical and chemi- 
cal evolution of a gas-rich dwarf galaxy. They found that most 
of the metals reside in the cold gas phase already after a few 
tens of Myr. This result is mainly due to the assumed low SNII 
heating efficiency, and justifies the generally adopted homoge- 
neous and instantaneous mixing of gas in chemical evolution 
models. On the other hand, Rieschick & Hensler (2000) pre- 
sented a chemodynamical model of the evolution of a dwarf 
galaxy in which most of the metals undergo a cycle lasting al- 
most 1 Gyr before being incorporated into the ISM. A sim- 
ilar result has been obtained by Tenorio-Tagle (1996). More 



recently, in the work of Harfst et al. (2006) a multi-phase de- 
scription of the ISM was implemented: the ISM was treated as 
two dynamically independent gas phases: a cloudy and a dif- 
fuse one. The two gas phases are linked by the processes of 
condensation/evaporation and a drag force. In this model, it is 
assumed that stars form in clouds and that clouds are destroyed 
by stellar feedback, thereby self-regulating the star formation. 
Each cloud is assumed to form stars only after a time of in- 
activity, because not all the gas in clouds is dense or unstable 
enough for immediate star formation. In this work the time of 
inactivity considered is a few hundred Myr. 

Here, we simply implemented our "metal cooling delay 
model" in the detailed chemical evolution model of F04 by 
means of a delay in the chemical enrichment. 

The paper is organized as follows: in Sect. 2 we describe 
the reference model of F04, in Sects. 3 and 4 we report the nu- 
cleosynthesis prescriptions and the observational data, respec- 
tively. In Sect. 5 we present our galactic fountains delay model, 
and in Sect. 6 the "metal cooling delay model". In Sect. 7 we 
report and discuss our results. Finally, we draw the main con- 
clusions in Sect. 8. 



2. The chemical evolution model for the Milky Way 

In the model of F04, the Galaxy is assumed to have formed by 
means of two main infall episodes: the first forms the halo and 
the thick disk, the second the thin disk. The time-scale for the 
formation of the halo-thick disk is 0.8 Gyr and the entire for- 
mation period does not last more than 2 Gyr. The time-scale 
for the thin disk is much longer, 7 Gyr in the solar vicinity, 
implying that the infalling gas forming the thin disk comes 
mainly from the intergalactic medium and not only from the 
halo (Chiappini et al. 1997). Moreover, the formation timescale 
of the thin disk is assumed to be a function of the galactocen- 
tric distance, leading to an inside out scenario for the Galaxy 
disk build-up (Matteucci & Francois 1989). The galactic thin 
disk is approximated by several independent rings, 2 kpc wide, 
without exchange of matter between them. 

The main characteristic of the two-infall model is an almost 
independent evolution between the halo and the thin disk (see 
also Pagel & Tautvaisienne 1995). A threshold gas density of 
7M Q pc~ 2 in the star formation process (Kennicutt 1989, 1998, 
Martin & Kennicutt 2001) is also adopted for the disk. 

The model well reproduces already an extended set of ob- 
servational constraints. Some of the most important observa- 
tional constraints are represented by the various relations be- 
tween the abundances of metals (C,N,0,a-elements, iron peak 
elements) as functions of the [Fe/H] abundance (see Chiappini 
et al. 2003a, b and F04) and by the G-dwarf metallicity distri- 
bution. It is worth mentioning here that, although this model 
is probably not unique, in this respect it reproduces the major- 
ity of the observed features of the Milky Way. Many of the 
assumptions of the model are shared by other authors (e.g. 
Prantzos & Boissier 2000, Alibes et al. 2001, Chang et al. 
1999). 
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The equation below describes the time evolution of G„ 
namely the mass fraction of the element ;' in the gas: 

Gi{r,t) = -Kr,tyXi{r,t) 



+ J iff(r, t - T m )Q mi (t - T m )(p(m)dm 



Mbm 

+A Ia J <f>(M B ) 



M Bm 



0.5 



J f(vMr,t-T m2 )Q™ Ia (t-T m2 )dLi 

[Mm 



dM„ 



Mbm 



+(1 - A Ia ) J t/f(r, t - T m )Q mi (t - T m )<f>{m)dm 

M Bm 

+ J t//(r, t - T m )Q mi (t - T m )4>(m)dm 

Mbm 



+ X A Mr,t) 



(1) 



where f) is the abundance by mass of the element ;' and 
Q m i indicates the fraction of mass restored by a star of mass 
m in form of the element i, the so-called "production matrix" 
as originally defined by Talbot & Arnett (1973). We indicate 
with Mi the lightest mass which contributes to the chemical 
enrichment and it is set at 0.8M o ; the upper mass limit, M v , is 
set at 1OOM . 

The star formation rate (SFR) if/(r, t) is defined as: 

>2(fc-l) ; 



if,(r,i) 



V (t-1 



Z(r,tGal) \ , . 

E(r,f) / gas 



(r,t), 



(2) 



where v is the efficiency of the star formation process and is set 
to be 2 Gyr 1 for the Galactic halo and 1 Gyr 1 for the disk. 
2(r, t) is the total surface mass density, S(r G , f) the total surface 
mass density at the solar position, G gas (r, f) the surface density 
normalized to the present time total surface mass density in the 
disk Z/)(r, t Ca i), where t Ga i = 14 Gyr is the age assumed for 
the Milky Way and r Q = 8 kpc the solar galactocentric dis- 
tance (Reid 1993). The gas surface exponent, k, is set equal to 
1 .5. With these values of the parameters, the observational con- 
straints, in particular in the solar vicinity, are well fitted. Below 
a critical threshold of the gas surface density (7M pc~ 2 ) we as- 
sume no star formation. This naturally produces a hiatus in the 
SFR between the halo-thick disk and the thin disk phases. 

For the IMF, we use that of Scalo (1986), constant in time 
and space. r m is the evolutionary lifetime of stars as a function 
of their mass m (Maeder & Maynet 1989). 

The Type la SN rate has been computed following Greggio 
& Renzini (1983) and Matteucci & Greggio (1986) and it is 
expressed as: 



M B . 



R 



SNela 



M Bm 



0.5 



= A Ia j <f>(M B ) j ' f(jl)lf/(t - T Ml )d^l 



dM R 



(3) 



where M 2 is the mass of the secondary, M B is the to- 
tal mass of the binary system, fi = M 2 /Mb, // m 
max [M 2 (t)/M B , (M B - 0.5M BM )/M B ], M Bm = 3M , M BM = 
16M Q . The IMF is represented by (p{M B ) and refers to the total 
mass of the binary system for the computation of the Type la 
SN rate, f(/i) is the distribution function for the mass fraction 
of the secondary: 



(4) 



with y = 2; Ai a is the fraction of systems in the appropriate 
mass range, which can give rise to Type la SN events. This 
quantity is fixed to 0.05 by reproducing the observed Type la 
SN rate at the present time (Mannucci et al. 2005). Note that 
in the case of the Type la SNe the"production matrix" is indi- 
cated with Q S J? la because of its different nucleosynthesis con- 
tribution (for details refer to Matteucci & Greggio 1986 and 
Matteucci 2001). 

The term A(r, t) represents the accretion term and is defined 



A(r, t) = a(r)e 



-t/T H (r) 



+ b{r)e 



-(t-W)/r D (r) 



(5) 



X Ai are the abundances in the infalling material, which is as- 
sumed to be primordial, while t max = 1 Gyr is the time for 
the maximum infall on the thin disk, t# = 0.8 Gyr is the time 
scale for the formation of the halo thick-disk and T D (r) is the 
timescale for the formation of the thin disk and is a function 
of the galactocentric distance (formation inside-out, Matteucci 
and Francois, 1989; Chiappini et al. 2001). 
In particular, we assume that: 



t d = 1.033r(kpc)- 1.267 Gyr 



(6) 



Finally, the coefficients a(f) and b(f) are obtained by imposing 
to fit the observed current total surface mass density in the thin 
disk as a function of galactocentric distance given by: 



cr(r) = cr D e- r/rD , 



(7) 



where cr D =531 M Q pc 2 is the central total surface mass den- 
sity and r D = 3.5 kpc is the scale lenght. 

3. Nucleosynthesis prescriptions 

For the nucleosynthesis prescriptions of the Fe and the other 
elements (namely O, S, Si, Ca, Mg, Sc, Ti, V, Cr, Zn, Cu, Ni, 
Co and Mn ), we adopted those suggested in F04. They com- 
pared theoretical predictions for the [el/Fe] vs. [Fe/H] trends 
in the solar neighborhood for the above mentioned elements 
and they selected the best sets of yields required to best fit the 
data. In particular for the yields of SNe II they found that the 
Woosley & Weaver (1995) (hereafter WW95) ones provide the 
best fit to the data. In fact, no modifications are required for 
the yields of Ca, Fe, Zn and Ni as computed for solar chemi- 
cal composition. For oxygen the best results are given by the 
WW95 yields computed as functions of the metallicity. For the 
other elements, variations in the predicted yields are required 
to best fit the data (see F04 for details). 

In particular, the Mg yields require some adjustments since 
Mg yields from massive stars are generally too low and predict 
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a a too low Mg solar abundance. However, it is worth noting 
that this fact does not affect the [Mg/Fe] vs. [Fe/H] diagram if 
the predicted ratios are normalized to the predicted solar abun- 
dances. What matters in these diagrams are the proportions 
of a given elements produced in different mass ranges: Mg is 
mainly produced in massive stars and therefore, irrespective of 
the predicted Mg abundance one always obtains an overabun- 
dance of Mg relative to Fe in metal poor stars followed by a 
decrease of the Mg/Fe ratio at the onset of Type la SNe which 
produce the bulk of Fe. In any case, F04 increased artificially 
the yields of Mg from WW95 to obtain also a good solar Mg 
abundance and to give an indication to nucleosynthesis model- 
ers about what might be necessary to fit the absolute Mg abun- 
dance. For the Fe-peak elements the situation is more complex 
and still very uncertain and in any case we do not show results 
for these elements in this paper. Concerning O, the best agree- 
ment with the [O/Fe] vs. [Fe/H] and the solar O abundance, as 
measured by Asplund et al. (2005), is obtained by adopting the 
original WW95 yields from massive stars as functions of metal- 
licity. The same is not true for Fe: it was already known since 
the paper of Timmes & al. (1995) that the Fe yields as func- 
tions of metallicity overestimate the solar Fe and many people 
use those yields divided by a factor of 2. Alternatively, one can 
use the yields for solar chemical composition for the whole 
galactic lifetime and the result is the same. This is due to the 
still existing uncertainties in the Fe yields. 

Concerning the yields from Type SNela, revisions in the 
theoretical yields by Iwamoto et al.(1999) are required for Mg, 
Ti, Sc, K, Co, Ni and Zn to best fit the data. The prescription 
for single low-intermediate mass stars are by van den Hoek 
& Groenewegen (1997), for the case of the mass loss param- 
eter which varies with metallicity (see Chiappini et al. 2003a, 
model5). 

4. Observational data 

We used the collection of data elaborated in the work of F04, 
who adopted a data sample for stars in the solar neighborhood 
spanning a metallicity range from -4 dex to solar. In particular, 
for the very metal poor stars ([Fe/H] between -4 and -3 dex), 
the results from the UVES Large Program " First Stars" (Cayrel 
et al. 2003) were adopted. 

For the abundances in the remaining range of [Fe/H], they 
adopted already published data in the literature from various 
sources: for all the elements studied in this paper (O, Mg, Si) 
F04 used the data of Carney et al. (1997), Nissen & Schuster 
(1997), Ryan et al. (1991), Edvardsson et al. (1993), Matteucci 
et al. (1993 and references therein). Moreover, concerning the 
Mg and Si, but not O, the collection of F04 includes the data 
of Stephens (1999), Carretta et al. (2002), Mc William et al. 
(1995), Fulbright (2000), Gratton & Sneden (1988). For the O, 
F04 also used the data of Nissen et al. (2002). Here, for the O, 
we also took into account the data by Bensby et al. (2004). 

All of these data are normalized to the solar abundances of 
Grevesse and Sauval (1998) with the exception of oxygen for 
which it is adopted the value of Asplund et al. (2005). 

For the galactic abundance gradient for the oxygen we used 
the set of data given by Andreievsky et al. (2002a-b) and Luck 



et al. (2003) by analyzing Galactic Cepheids (see Cescutti at al. 
2007). 

5. The galactic fountain delay model 

In Paper I we used the Kompaneets (1960) approximation for 
the evolution of the superbubble driven by sequential super- 
nova explosions. We assumed that supershells are fragmented 
by Rayleigh-Taylor instabilities, and considered each fragment 
as a cloud with an initial velocity given by the supershell ve- 
locity at the moment of fragmentation. Then, we followed the 
orbit of the clouds either ballistically (the cloud is subject only 
to the gravitational potential of the Galaxy) or taking into ac- 
count viscous interaction between the extra-planar gas halo and 
the cloud. 

In Tables 2-4 reported in Paper I, are summarized the re- 
sults for fragmentation times (f//, M ;) and the velocities (v„) of 
the supershells in the direction perpendicular to the galactic 
plane at those times for initial throwing coordinates fixed at 
4 kpc, 8 kpc and 12 kpc. In the Table Q] here we report the 
average time < t, ota i > necessary to create a cloud from a su- 
pershell plus the time necessary to the cloud to fall back into 
the Galactic disk: 

< hotal >= tfinal+ < t orBit >, (8) 

where < t or ti, > is the average time required for a cloud to 
return to the galactic disk once it leaves the supershell. These 
average values are calculated in the case of OB associations 
with 10, 50, 100, 500 SNe at 3 galactocentric distances: 4 kpc, 
8 kpc, 12 kpc. In Table [2] we show the maximum delay time 
hotal , at a fixed initial galactocentric distance of ejection. 

Table 1. The average time < t tota i > [Myr]. 





4 kpc 


8 kpc 


12 kpc 


10 SNe 


43 


53 


75 


50 SNe 


36 


54 


87 


100 SNe 


36 


57 


96 


500 SNe 


38 


75 


133 



Table 2. The maximum t ma i as a function of the galactocentric 
distance. 



4 kpc 


8 kpc 


12 kpc 


Maximum delay 48 Myr 


1 14 Myr 


245 Myr 



As mentioned in Sect.l, a Galactic fountain can affect the 
chemical evolution of the Galaxy only by means of a time- 
delay due to the non-negligible time taken by fountains to orbit 
around and fall back into the Galaxy. We implemented this de- 
lay in the chemical evolution model of F04 only during the thin 
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O 20 - 



m 10 



"1 1 <~ 

4 kpc 
8 kpc 
12 kpc 




t [Gyr] 

Fig. 1. The star formation rate expressed in M Q pc -2 Gyr -1 as 
a function of the galactic time (in Gyr), as predicted by the 
two infall model of Francois et al. (2004) at 4 (dotted line), 8 
(dashed line) and 12 kpc (solid line). 



Table 3. Beginning of the star formation in the disk after the 
halo-thick disk transition (?sfr)- 



halo-thick disk transition. In Fig. Q] we show the predicted star 
formation rate (SFR) expressed in M Q pc~ 2 Gyr -1 as a func- 
tion of the galactic time (Gyr) using the model of F04. We note 
that the SFR is the same for all galactocentric distances during 
the halo phase. The effect of the two infall model in presence 
of the threshold in the star formation is clearly shown: the star 
formation halts during the halo-thin disk transition and the du- 
ration of the gap varies as a function of galactocentric distance. 
At the solar neighborhood the gap is ^0.3 Gyr. This gap seems 
to be observed in the abundance ratio distribution as shown in 
Gratton et al. (2000) and Fuhrmann (1999). In Table[3]we re- 
port the times at different galactocentric distances at which the 
SFR in the thin disk begins (tsFu) in the F04 model. Ejecta 
of intermediate-mass stars and Type la SNe (which are most 
probably not clumped in OB associations) are not subject to 
the galactic fountain delay. 

As shown in Table [2] a single superbubble can produce a 
delay in the mixing up to ~ 250 Myr, but the average values are 
smaller, spanning the range between 40 - 140 Myr (Table [TJ. 
In the simulation we considered different values up to 1 Gyr. It 
is possible that the Milky Way formed super-star clusters dur- 
ing more intense phases of star formation, as observed in M82 
(Melo et al. 2005). A delay of 1 Gyr can then be obtained if the 
number of massive stars in a cluster is ~ 10 4 , which is consis- 
tent with the mass of the largest super-star clusters observed in 
M82. In any case, we considered this long delay as an extreme 
case. 



4 kpc 


8 kpc 


12 kpc 


Beginning of SF in the disk 1 . 1 Gyr 


1.3 Gyr 


1.7 Gyr 



disk phase. In fact, the dynamical effect of a galactic fountain 
requires the presence of a stratified ISM, which arises only af- 
ter the end of the halo phase. In the past, several authors (i.e. 
Mac Low et al. 1989 and Tenorio-Tagle 1996), studied the su- 
pershell evolution produced by a sequential explosion of Type 
II SNe. They emphasized that a break-out event, necessary for 
a galactic fountain, requires that the OB associations must sit 
on a plane stratified disk where the density decreases along the 
Z direction. In the work of Mac Low et al. (1989) a scale height 
of 100 pc for the ISM was considered. This is the reason why 
OB associations existed at early times during the formation of 
the halo but did not produce galactic fountains. Galactic foun- 
tains could have formed perhaps during the thick disk phase. 
However, by computing a fountain with a scale height of 1 kpc 
(typical of the thick disk), the fountains start at 3 kpc, i.e. 
well within the halo. At this height the stratification is much 
weaker and the ISM temperature is high, therefore the sound 
speed is high and the shock can be easily dissipated through 
sound waves. Moreover the thick disk phase lasts for a very 
short time (see Fig. [TJ. We can then conclude that (hypotheti- 
cal) thick disk fountains would have a negligible impact on the 
chemical evolution of the Galaxy. 

Therefore, we consider the galactic fountain delay only for 
ejecta from massive stars (with mass > 8 M ) born after the 



6. The metal cooling delay model 

Roy & Kunth (1995) studied the dispersal and mixing of oxy- 
gen in the ISM of gas-rich galaxies. Due to a variety of hy- 
drodynamical processes operating at scales ranging from 1 pc 
to greater than 10 kpc in the ISM, there are difficulties in esti- 
mating the dispersal and mixing timescale: (/) on large galactic 
scales (1 > / > 10 kpc), turbulent diffusion of interstellar 
clouds in the shear flow of galactic differential rotation is able 
to wipe out azimuthal O/H fluctuations in less than 10 9 yrs; (//) 
at the intermediate scale (100 > I > 1000 pc), cloud collisions 
and expanding supershells driven by evolving associations of 
massive stars, differential rotation and triggered star formation 
will re-distribute and mix gas efficiently in about 10 8 yrs; (Hi) 
at small scales (1 > I > 100 pc), turbulent diffusion may be 
the dominant mechanism in cold clouds, while Rayleigh-Taylor 
and Kelvin-Helmhotz instabilities quickly develop in regions of 
gas ionized by massive stars, leading to full mixing in < 2 x 10 6 
yrs. Malinie et al. (1993) claimed that due to chemical inhomo- 
geneities observed for example in the age-metallicity relation 
in the solar neighborhood, re-mixing and star formation may 
be delayed by 10 8 ~ 9 yr. 

Following the work of Thomas et al. (1998) we relaxed the 
IMA in the chemical evolution model by splitting the gas com- 
ponent into two different phases (cold and warm). The gas in 
the interstellar medium is heated by SN events and the stars 
cannot form in this "warm gas phase"; only after the gas cools 
stars can form and then pollute the ISM. We implemented the 
two gas phases simply by means of a delay in the chemical en- 
richment. In this case: i) all stellar masses (both Type II and 
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Type la SNe) and ii) both halo and disk stars, are affected by 
this delay. 

7. Results 

In this section we discuss our results. First of all we show how a 
galactic fountain delay can affect the chemical evolution model, 
then we present the results for the metal cooling delay model. 

We show the model predictions for O, Mg, Si and Fe ob- 
tained in particular for the relations [el/Fe] versus [Fe/H] com- 
pared with the observational data and try to put constraints on 
the maximum possible delay. 

7.1. Results for the galactic fountain delay model 

We have taken into account different values for the delay orig- 
inated by galactic fountain events from sequential SN explo- 
sions: 0.1 Gyr, 0.2 Gyr, 0.5 Gyr, 1 Gyr. We tested the effect at 
different galactocentric distances: 4, 8, 12 kpc. 

Concerning the solar neighborhood, in Table [4] we report 
the solar abundances by mass as predicted by the F04 model 
and the delay models compared to the solar values of Asplund 
et al. (2005). In Fig. [2] we show the results for [O/Fe] versus 
[Fe/H] in the solar neighborhood compared to the observational 
data described in Sect 4. Our results are normalized to the "best 
model" of F04 values at 9.5 Gyr (the age of the formation of the 
solar system). As we can see from this figure and from Table|4] 
the solar abundances predicted by our models are substantially 
in agreement with those derived by Asplund et al. (2005). The 
standard model with no delay shows a discontinuity at ~ -0.8 
dex, and this is due to the gap in the SFR discussed before. In 
fact, during the gap O stops being produced whereas Fe keeps 
being produced by Type la SNe. The effect of the delay due to 
the fountains is therefore to increase this effect by increasing 
the duration of the gap. 

In the description of the delay model we emphasized that 
galactic fountains can be seen only in disk stars. The typical 
delay, due to the fact that massive stars are clumped in OB as- 
sociations, is 0.1 Gyr. Referring to Fig. [2] we present a first 
strong result: a delay of 0.1 Gyr produces a negligible effect 
on the chemical evolution of the Galaxy in the solar neighbor- 
hood. We note that the data show a spread in the halo-thick disk 
transition phase ([Fe/H]^-0.77). This spread can be explained 
in terms of a combination of: i) a threshold in the star forma- 
tion, ii) a delay in the chemical enrichment from the massive 
stars due to galactic fountain effect. Moreover, we note that the 
maximum possible delay, in order not to break the agreement 
with data, must be lower than 1 Gyr. 

The same effect is shown in Fig. [3] and in Fig.|4]where we 
analyzed the behavior of such a delay in a chemical evolution 
model for other two a elements: Mg and Si. 

It is also interesting to test the fountain delay at different 
galactocentric distances. In Figs. [5] [6] and [7] we report our re- 
sults for O, Mg and Si, respectively. In each figure we plot the 
[el/Fe] vs. [Fe/H] at 4 and 12 kpc galactocentric distances. We 
note that the results are strongly dependent on the galactocen- 
tric distance. In fact, as explained before, the drop of the [el/Fe] 
quantity in the [el/Fe] vs. [Fe/H] is due to the iron ejected by 




Fig. 2. The abundance ratio [O/Fe] as a function of [Fe/H] in 
the solar neighborhood. We compared the Francois et al. (2004) 
model (solid line) with our galactic fountain delay models: dot- 
ted line (0. 1 Gyr), short dashed line (0.2 Gyr), long dashed line 
(0.5 Gyr), dashed dotted line (1 Gyr). The data are taken from 
the collection used in Francois et al. (2004) (filled green trian- 
gles) and from Bensby et al. (2005) (filled blue circles). 



8 kpc 



0.5 



-0.5 - 




no delay 

0.1 Gyr 

0.2 Gyr . 
__ 0.5 Gyr 
1 Gyr 



-4 -3 -2 -1 

[Fe/H] 

Fig. 3. The abundance ratio [Mg/Fe] as a function of [Fe/H] in 
the solar neighborhood. Notation as in Fig. [2] 



the Type la SNe. The discontinuity predicted at 4 kpc is larger 
than that predicted at larger galactocentric distances. This ef- 
fect is due to the higher specific Type la SN rate at 4 kpc than 
at larger distances. The specific Type la SN rate is the SNIa rate 
per unit mass of gas and it increases toward the Galactic center, 
the reason being that while the Type la SN rate is a factor of 
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Table 4. Mass fractions of several elements yielded by the Francois et al. (2004) model and the delay model at the galactic age 
of 9.5 Gyr compared to the solar values of Asplund et al. (2005). 











Asplund et al. (2005) 


Francpis et al. (2004) 


0.1 Gyr 


0.2 Gyr 


0.5 Gyr 


1 Gyr 


He 


2.55 xlO- 1 


2.52 xlO-' 


2.52 xlO- 1 


2.52 xlO" 1 


2.52 xlO-' 


2.52 xlO- 1 


C 


2.21 xl0~ 3 


2.03 xl0~ 3 


2.05 xl0~ 3 


2.06 xl0~ 3 


2.11 xlO" 3 


2.17 xlO" 3 


N 


6.32 xl0~ 4 


9.73 xl0~ 4 


9.75 xlO" 4 


9.77 xl0~ 4 


9.81 xlO- 4 


9.85 xlO" 4 


O 


5.48 xl0~ 3 


5.27 xl0~ 3 


5.31 xlO" 3 


5.35 xlO" 3 


5.46 xlO" 3 


5.62 xlO" 3 


Mg 


6.17 xl0~ 4 


6.66 xl0~ 4 


6.70 xl0~ 4 


6.74 xlO" 4 


6.86 xlO 4 


7.03 xlO" 4 


Si 


6.80xl0~ 4 


8.14 xl0~ 4 


8.17 xlO" 4 


8.20 xlO" 4 


8.30 xlO" 4 


8.43 xlO" 4 


S 


3.31 xl0~ 4 


3.95xl0- 4 


3.97xl0- 4 


3.98xl0 4 


4.03 xlO 4 


4.09xl0- 4 


Ca 


6.13 xl0~ 5 


5.40 xl0~ 5 


5.42 xlO" 5 


5.44 xlO- 5 


5.50xlO- 5 


5.58 xlO- 5 


Fe 


1.18 xl0~ 3 


1.14 xl0~ 3 


1.15 xlO" 3 


1.15 xlO" 3 


1.15 xlO" 3 


1.16 xlO" 3 




8 higher than that in the solar vicinity, because of the higher 
SFR, the mass of gas is higher by only a factor of 2 relative to 
the solar region. Therefore, the effect of the pollution by Type 
la SNe during the gap is enhanced. 

Again in Figs.[5]|6]and[7]we note that, at a fixed galactocen- 
tric distance, the effect of the galactic fountain delay depends 
also on the considered element: the Si which is also produced 
by Type la SNe in a non negligible amount shows a smaller 
drop of the [Si/Fe] quantity compared to O and Mg. 

Even if the galactic fountain effect in the [el/Fe] values is 
larger at 4 kpc because of the higher rate of Type la SNe, the 
delays expected at 4 kpc are smaller with respect to outward 
galactocentric distances (as shown Tabs Q] and [2]), as a conse- 
quence of the higher gravitational pull. In Fig. [9] we show the 
[O/Fe] vs [Fe/H] at 4, 8 and 12 kpc considering the maximum 
delays in each radius. In the case of a delay produced by a sin- 
gle OB association, the effect of galactic fountains is negligible. 



Finally we explored the effects of galactic fountains on the 
abundance gradient for the oxygen. In Fig.[10]we plot the [O/H] 
as a function of the galactocentric distance. We note that the av- 
erage delay due to galactic fountains is longer than 0. 1 Gyr only 
for a OB association composed of 500 SNe at 12 kpc (Table[TJ; 
this event is very unlikely, then the delay of 0.1 Gyr can be 
considered like a average upper limit for a delay produced by a 
single OB association at a galactocentric distance in the range 
4-12 kpc. Referring to Fig[l0]we can conclude that the time 
delay produced by a Galactic fountain originated by an OB as- 
sociation has a negligible effect on the abundance gradient in 
the Galaxy disk. 

7.2. Results for the metal cooling delay model 

In this section we present the results for the metal cooling delay 
model. Referring to Malinie et al. (1993) we considered mixing 
delays of 0.5 Gyr and 1 Gyr. In Figs.[TT|[T2]and[T3]we snow 
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-0.2 



no delay 

0.1 Gyr 

0.2 Gyr 

0.5 Gyr 

1 Gyr 



-2 -1 
[Fe/H] 



Fig. 6. The abundance ratio [Mg/Fe] as a function of [Fe/H] at 
4 and 12 kpc. Notation as in Fig. [2] 



o.oi 



4 kpc 
8 kpc 
12 kpc 



10 



15 



t [Gyr] 



Fig. 8. The Type la SN rate expressed in pc~ 2 Gyr 1 as a func- 
tion of the galactic time (Gyr), as predicted by the two infall 
model of Francois et al. (2004) at 4, 8 and 12 kpc. 




0.6 
^0.4 
^0.2 



no delay 

0.1 Gyr 

0.2 Gyr 

_ - 0.5 Gyr 

- 1 Gyr 



-0.2 



-2 -1 
[Fe/H] 



Fig. 7. The abundance ratio [Si/Fe] as a function of [Fe/H] at 4 
and 12 kpc Notation as in Fig. [2] 




Fig. 9. The abundance ratio [O/Fe] as a function of [Fe/H] at 4, 
8 and 12 kpc for relative maximum delays as reported in Tab|2] 



the effect of the metal cooling predicted by our model using 
the [el/Fe] vs. [Fe/H] plots. The work of Thomas et al. (1998) 
already discussed the effect of a two- phase gas model relative 
to Mg. Referring to Fig. [12] we obtain a similar result, but the 
effect of our delay in the mixing is smaller. Also for the Si we 
have a similar behavior as can be seen in the Fig. [13] The dif- 
ferences between ours and Thomas et al.'s results are probably 
due to the different adopted yields for massive stars and IMF. 
In fact, Thomas et al. adopted Thielemann et al. (1996) yields 
and the Salpeter (1955) IMF. 



A different result is obtained for the [O/Fe] vs [Fe/H] re- 
ported in Fig. QT] In this case, we find that the delay in the 
mixing leads to big differences in terms of chemical evolution 
models. We see that in the halo-thick disk phase, the predicted 
[O/Fe] values are much smaller than for the model with the 
IMA. The reason for this resides in the fact that in the model of 
F04 the yields used for the O are the WW95 metallicity depen- 
dent ones and in this case the delay induces a situation where 
the yields for Z=0, lower than those computed for the other 
metallicities, act for a longer time, thus producing the large de- 
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Fig. 10. Predicted and observed O abundance gradients in the 
galactocentric distance range 4-14 kpc. The data points are 
from Cepheids. The large squares with error bars represent av- 
erage of the points with their errors. The black solid line is the 
Francois et al. (2004) model normalized to the mean value of 
the bin centered at 8 kpc (see Cescutti et al. 2007). 




-4 -3 -2 -1 

[Fe/H] 



Fig. 11. The abundance ratio [O/Fe] as a function of [Fe/H] 
in the solar neighborhood predicted by our metal cooling de- 
lay model. For the O we considered the metallicity dependent 
yields of WW95. Solid line: model without delay; dotted line 
and dashed dotted line: 0.5 and 1 Gyr delays, respectively. 
Symbols as in Fig. [2] 



pression in the [O/Fe] ratio. To demonstrate that this result is 
solely due to the choice of nucleosynthetic yields, we plot in 
Fig. [14] the metal cooling delay model results using, for the 
oxygen, the yields at solar metallicity given by WW95. As we 
can see from this figure, in this case the models with delay do 
not differ considerably compared to the reference model and 
the oxygen behaves similarly to the other c-elements. In Fig. 
[TBI we report the [el/Fe] vs. [Fe/H] relations for the Mg and Si 
using for these two elements the metallicity dependent yields of 
WW95. As expected, the [el/Fe] ratios show a large depression 
in the halo-thick disk phase as in Fig. [TTIdue to the choice of 
metallicity dependent yields. Therefore we can conclude that 
the metallicity-dependent yields of WW95 are not compatible 
with a delayed enrichment of the halo. 

Finally, in Figs. [16] and [17] we report age-metallicity rela- 
tions in the solar neighborhood for [Fe/H] and [Mg/H] respec- 
tively, predicted by our metal cooling delay model using the 
stellar yields adopted in F04. Concerning the evolution of the 
[Fe/H] we compare our results with the observational data of 
Ramirez et al. (2007). As one can see, the effects of the de- 
lays are more visible in the age-metallicity relations than in the 
abundance ratios versus [Fe/H], and the reason is that in the 
abundance ratios these differences are partly erased. We con- 
clude that all delays considered are compatible with the obser- 
vations owing to the very large data spread for the age metal- 
licity relation, but also that the age-metallicity relation is not a 
good constraint for chemical models. 



0.5 



-0.5 



8 kpc 




no delay 
0.5 Gyr 
1 Gyr 



-4 -3 -2 -1 

[Fe/H] 

Fig. 12. The abundance ratio [Mg/Fe] as a function of [Fe/H] 
in the solar neighborhood predicted by our metal cooling delay 
model. Notation as in FigJTTI 



8. Conclusions 

In this paper we have studied the relaxation of IMA in chem- 
ical evolution models by means of the effects of a delay in 
the chemical enrichment produced by galactic fountain events 
and by means of a delay in the mixing due to chemical inho- 
mogeneities in the disk (our metal cooling delay model). The 
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Fig. 13. The abundance ratio [Si/Fe] as a function of [Fe/H] 
in the solar neighborhood predicted by our the metal cooling 
delay model. Notation as in FigfTTI 



Fig. 15. The abundance ratio [el/Fe] as a function of [Fe/H] 
for the Si and Mg in the solar neighborhood predicted by our 
metal cooling delay model. For the Si and Mg we considered 
the metal dependent yields of WW95. Notation as in FigfTTI 



r-0.5 




0.5 



no delay 

0.5 Gyr 

_.. 1 Gyr 



-4 -3 -2 -1 

[Fe/H] 

Fig. 14. The abundance ratio [O/Fe] as a function of [Fe/H] in 
the solar neighborhood predicted by our metal cooling delay 
model. For the O we considered the solar yields of WW95. 
Notation as in FigfTTI 

main purpose of this work was testing the results of Paper I 
in a chemical evolution model. In the present paper we tested 
the "delay" effect due to the finite and not negligible time for a 
cloud spend to fall back onto the disk. 

Our main conclusions can be summarized as follows: 

- In the solar neighborhood we showed that the average delay 
produced by the galactic fountains has a negligible effect on 
the chemical evolution for all the a elements we studied. 



no delay 
0.5 Gyr 

1 Gyr 




15 



t [Gyr] 



Fig. 16. The abundance ratio [Fe/H] as a function of the galactic 
time [Gyr] in the solar neighborhood predicted by our metal 
cooling delay model. The data are taken from Ramirez et al. 
(2007) (filled red triangles). Notation as in FigfTTI 



- In the [el/Fe] versus [Fe/H] relations, the main feature of 
the galactic fountain is an enhancement of the drop in the 
[el/Fe] ratios occurring because of the two infall scenario. 
In fact, the drop in the standard model of F04 is due to the 
halt of the SFR which produces a halt in the production of 
the Q'-elements while the Fe from Type la SNe continues 
to be produced. The galactic fountain delay has the effect 
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Fig. 17. The abundance ratio [Mg/H] as a function of the galac- 
tic time [Gyr] in the solar neighborhood predicted by our metal 
cooling delay model. Notation as in FigfTTI 

of increasing the period during which there is no pollution 
from Type II SNe. 

- Results produced by the model with a galactic fountain de- 
lay of 1 Gyr are not compatible with observational data. 
On the other hand, any delay < 1 Gyr could be acceptable, 
given the observed spread in the data. 

- The time delay produced by a Galactic fountain originated 
by a OB association has a negligible effect on the abun- 
dance gradients in the Galaxy disk. 

- The metal cooling delay model with the assumed delays 
(which are those suggested in the literature), has a very 
small effect on the chemical evolution in the solar neigh- 
borhood if yields not depending on metallicity are used. 

- On the other hand, in the case of the metal dependent yields 
of WW95, the results differ substantially from the reference 
model of F04 and do not fit the observations. 
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